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ABSTRACT 


The  assumption  of  self-preservation  allows  for  an  analytical  determination  of  the 
energy  decay  in  isotropic  turbulence.  Batchelor  (1948),  who  was  the  first  to  carry  out  a 
detailed  study  of  this  problem,  based  his  analysis  on  the  assumption  that  the  Loitsianskii 
integral  is  a  dynamic  invariant  —  a  widely  accepted  hypothesis  that  was  later  discovered 
to  be  invalid.  Nonetheless,  it  appears  that  the  seif-preserving  isotropic  decay  problem 
has  never  been  reinvestigated  in  depth  subsequent  to  this  earlier  work.  In  the  present 
paper  such  an  analysis  is  carried  out,  yielding  a  much  more  complete  picture  of  self¬ 
preserving  isotropic  turbulence.  It  is  proven  rigorously  that  complete  self-preserving 
isotropic  turbulence  admits  two  general  types  of  asymptotic  solutions:  one  where  the 
turbulent  kinetic  energy  K  ~  and  one  where  K  ~  t~°‘  with  an  exponent  a  >  1  that  is 
determined  explicitly  by  the  initial  conditions.  By  a  fixed  point  analysis  and  numerical 
integration  of  the  exact  one-point  equations,  it  is  demonstrated  that  the  K  ~  power 
law  decay  is  the  asymptotically  consistent  high-Reynolds-number  solution;  the  K  ~ 
decay  law  is  only  achieved  in  the  limit  as  f  oo  and  the  turbulence  Reynolds  number 
vanishes.  Arguments  are  provided  which  indicate  that  a  K  ~  i~^  power  law  decay  is  the 
asymptotic  state  toward  which  a  complete  self-preserving  isotropic  turbulence  is  driven  at 
high  Reynolds  numbers  in  order  to  resolve  the  imbalance  between  vortex  stretching  and 
viscous  diffusion.  Unlike  in  previous  studies,  the  asymptotic  approach  to  a  self-preserving 
state  is  investigated  which  uncovers  some  surprising  results. 


•Research  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA  Contract 
No.  NASl-1860.5  while  the  author  was  in  re.sidence  at  the  Institute  for  Computer  Applications  in  Science 
and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  V\  ‘23G6.')-.u22.o. 
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1.  INTRODUCTION 


Despite  the  fact  that  isotropic  turbulence  constitutes  the  simplest  type  of  turbulent 
flow,  it  is  still  not  possible  to  render  the  problem  analytically  tractable  without  the 
introduction  of  additional  hypotheses.  The  idealization  of  self-preservation  —  wherein 
the  two-point  double  and  triple  longitudinal  velocity  correlations  are  assumed  to  admit 
self-similar  solutions  with  respect  to  a  single  length  scale  L{t)  —  has  served  as  a  useful 
hypothesis  since  its  introduction  by  von  Karman  and  Howarth  (1938).  In  another  classic 
paper  that  followed,  Batchelor  (1948)  studied  the  energy  decay  in  self-preserving  isotropic 
turbulence  in  considerable  detail.  He  concluded  that  the  only  complete  self- preserving  so¬ 
lutions  that  were  internally  consistent  existed  at  low  turbulence  Reynolds  numbers  where 
the  turbulent  kinetic  energy  K  ~  —  a  power  law  consistent  with  the  final  period 

of  decay.  Batchelor  (1948)  also  found  a  self-preserving  solution  to  the  Karman-Howarth 
equation  in  the  limit  of  infinite  Reynolds  numbers  for  which  Loitsianskii’s  integral  was  an 
invariant.  This  solution  —  wherein  L{t)  is  the  integral  length  scale  A  and  K  ~  — 

was  put  forth  by  Batchelor  as  the  only  complete  self-preserving  solution  at  high  Reynolds 
numbers.  Of  course,  additional  partial  self-preserving  solutions  were  shown  by  Batchelor 
to  exist  in  other  Reynolds  number  regimes. 

Objections  were  later  raised  against  the  use  of  the  Loitsianskii  integral  as  a  dynamic 
invariant;  at  high  Reynolds  numbers  this  integral  can  be  shown  to  be  a  weak  function  of 
time  (see  Proudman  and  Reid  1954  and  Batchelor  and  Proudman  1956).  Saffman  (1967) 
proposed  an  alternative  dynamic  invariant  which  yields  a  K  ~  power  law  decay  in 
the  infinite  Reynolds  number  limit  (see  Hinze  1975).  While  the  results  of  Batchelor  and 
Saffman  formally  constitute  complete  self-preserving  solutions  to  the  inviscid  Karman- 
Howarth  equation,  it  must  be  kept  in  mind  that  they  only  exhibit  partial  self-preservation 
with  respect  to  the  full  viscous  equation.  Namely,  there  is  self-preservation  only  for  the 
range  of  energy  containing  eddies  with  integral  length  scale  A  (here,  A/A  <<  1  where  A  is 
the  Taylor  microscale).  These  two  solutions  have  been  widely  accepted  in  the  turbulence 
literature  as  the  predicted  decaj  laws  for  sclf-pieseiving  isotropic  turbulence  at  high 
Reynolds  numbers. 

Implicit  in  the  analysis  of  Batchelor  (1948)  is  the  existence  of  a  complete  self-preserving 
solution  consistent  with  high  Reynolds  numbers  —  namely  a  K  ~  power  law  decay. 
The  collapsing  length  scale  L{t)  for  this  full  self-preserving  solution  is  necessarily  the 
Taylor  microscale  (i.e.,  for  any  complete  self-preserving  solution  of  the  viscous  Karman- 
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Howarth  equation  we  must  have  L  oc  A).  This  solution  ~  which  was  postulated  earlier  by 
Dry  den  (1943)  —  was  dismissed  by  Batchelor  on  the  grounds  that  Loitsianskii’s  integral 
was  not  a  dynamic  invariant.  While  this  result  by  Dry  den  has  been  mentioned  subse¬ 
quently  in  the  literature  (c.f.  Hinze  1975,  Monin  and  Yaglom  1975  and  Korneyev  and 
Sedov  1976),  it  has  largely  been  disregarded  by  tlie  turbulence  community.  The  reason 
for  this  appears  to  be  two-fold:  (a)  a  K  ~  power  law  decay  has  not  been  observed 
in  the  most  accurate  isotropic  decay  experiments,  and  (b)  since  A/A  — >  0  as  /{e  — >  oo, 
doubts  can  be  cast  on  the  suitability  of  the  Taylor  microscale  as  the  collapsing  length 
scale  of  the  energy  containing  eddies. 

Recently,  George  (1987,  1989)  revived  this  issue  concerning  the  existence  of  complete 
self-preserving  solutions  in  isotropic  turbulence.  In  an  interesting  paper  he  claimed  to 
find  a  complete  self-preserving  solution,  valid  for  all  Reynolds  numbers,  in  which  the 
kinetic  energy  decayed  as  K  ~  t”"  with  a  determined  by  the  initial  conditions.  George 
(1987)  —  who  based  his  analysis  on  the  dynamic  equation  for  the  energy  spectrum  rather 
than  on  the  Karman-Howarth  equation  —  made  no  explicit  mention  of  the  complete  self¬ 
preserving  K  ~  solution.  Strictly  speaking,  the  solution  presented  by  George  was  an 
alternative  seif-preserving  solution  to  that  of  Karman  and  Howarth  (1938)  and  Batchelor 
(1948)  since  he  relaxed  the  constraint  that  the  triple  longitudinal  velocity  correlation  be 
self-similar  in  the  classical  sense. 

The  purpose  of  the  present  paper  is  to  address  the  issue  of  complete  self-preservation 
in  an  effort  to  clarify  the  following  basic  questions: 

(1)  What  is  the  complete  self-preserving  solution  for  isotropic  turbulence  at  high 
Reynolds  numbers, 

(2)  What  detailed  predictions  does  this  solution  yield  for  the  energy  decay,  particularly 
during  the  initial  approach  to  a  state  of  complete  self-preservation,  and 

(3)  Is  this  solution  compatible  with  the  results  of  physical  experiments  and  alternative 
theoretical  approaches? 

In  so  far  as  tlie  first  two  points  are  concerned,  it  will  be  shown  unequivocally  that  the 
onl}.  coniplete  self-preserving  solution  that  is  consistent  with  a  high-Reynolds-number 
isotropic  turbulence  has  a  A  ~  t  '  asymptotic  power  law  decay.  Unlike  previous  studies, 
this  is  demonstrated  in  a  straightforward  manner  based  on  a  fixed  point  analysis  of  the 
one-point  equations.  This  analysis  leads  to  an  interesting  interpretation  of  the  physical 
significance  of  a  A  ~  t  ‘  powei  law  decay  aim  auuvva  us  to  examine  smail  departures 
from  a  state  of  complete  self-preservation.  The  detailed  jjredictions  of  this  complete  self- 
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preserving  solution  —  which,  to  the  best  of  our  knowledge  have  never  been  examined  in 
depth  in  the  literature  —  will  be  compared  with  the  results  of  physical  experiments  and 
alternative  theoretical  approaches  in  the  sections  to  follow. 


2.  THEORETICAL  BACKGROUND 


We  will  consider  isotropic  turbulence  governed  by  the  incompressible  Navier-Stokes 
equations 

r}ii  fin 

(1) 


dui  du,  _  dp  ^ 


OXi 


(2) 


where  Ui  is  the  velocity  vector,  p  is  the  pressure  and  i>  is  the  kinematic  viscosity.  The  two 
point  double  and  triple  longitudinal  velocity  correlations,  denoted  by  f{r,t)  and  k{r,t), 
respectively,  are  defined  in  the  standard  way: 


u(x,  <)«(x  +  r,t) 


k{r,t)  = 


ii2(x,f)n(x  +  r,t) 


(3) 


(4) 


(n2)3/2 

where  u  is  any  component  of  the  velocity,  x  and  x  +  r  are  aiiy  two  spatial  points  separated 
by  a  distance  r  =^j  r  j  in  the  direction  of  u,  and  an  overbar  denotes  a  spatial  average  (c.f. 
Hinze  1975).  For  isotropic  turbulence,  /  and  k  satisfy  the  Karman-Howarth  equation 

d{u^ 


-(321 

^,.2  J.  Qj, 


(5) 


which  is  obtained  directly  from  the  Navier-Stokes  equations.  The  turbulent  kinetic  energy 
K  =  ^UiUi  is  a  solution  of  the  differential  equation 


K  =  -< 


where 


_  ..  .  .2 


(6) 


(7) 


c  =  r  u)iUi  =  V  u: 

is  ine  turbulent  dissipaton  rate,  cu,  is  the  vorticity  vector  and  is  the  enstrophy.  The 
lurbulc'’*^  n'udi  rate  a  solution  of  the  dilfcientiai  equation 

7  „  7 


(  ! — t  (  e 

e  =  -  7-6’- 

3^15  ^  A  15  A 


(8) 
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where 


are  the  velocity  derivative  skewness,  the  coefficient  of  the  enstrophy  destruction,  the 
turbulence  Reynolds  number,  and  the  Taylor  microscale,  respectively.  Equations  (6)  and 
(8)  —  which  are  rearranged  versions  of  those  that  appear  in  Karman  and  Howarth  (1938) 
and  Batchelor  (1948)  —  are  obtained,  respectively,  by  setting  r  =  0  in  (5)  and  in  the 
equation  obtained  by  differentiating  (5)  twice  with  respect  to  r.  Since  Sk  =  SK{t)  and 
G  =  G{t)  are  directly  related  to  the  correlations  /  and  k  (which  cannot  both  be  obtained 
from  the  Karman- Howarth  equation  (5))  it  is  clear  that  the  problem  of  isotropic  decay 
is  not  closed.  In  order  to  achieve  closure,  additional  assumptions  must  be  made  such  as 
the  one  of  complete  self-preservation  that  will  be  discussed  in  this  paper. 

For  an  isotropic  turbulence  to  be  self-preserving  in  the  sense  of  Karman  and  Howarth 
(1938)  and  Batchelor  (1948),  we  must  have 


/(r,()  =  /(r/L)  (12) 

k(r,t)  =  k{r/L)  (13) 

where  L  =  L{t)  is  a  uniquely  specified  similarity  length  scale.  For  it  to  exhibit  complete 
self-preservation,  all  scales  of  the  turbulence  —  namely,  the  full  range  of  0  <  r  <  oo 
—  must  decay  according  to  (12)  and  (13)  (partial  self-preservation  is  satisfied  if  (12)  - 
(13)  only  apply  to  some  restricted  range  of  0  <  r  <  rmax)-  We  will  focus  our  attention 
on  complete  self-preserving  solutions  in  the  analysis  to  follow.  In  view  of  the  identity 
(Batchelor  1948) 

£  =  -lOuK 


it  follows  from  (11)  that 


d'^j 


dr'i 


8V1 

^  J  r=0 

(14) 

=  -1. 

(15) 

r=0 
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(16) 


Hence,  for  any  complete  self-preserving  isotropic  turbulence,  we  must  have 

^/"(O)  =  -1 

from  which  it  can  be  concluded  that 

L<x\  (17) 

since  f"{^)  is  a  constant.  It  therefore  follows  that  the  Taylor  microscale  is  the  only 
similarity  length  scale  that  can  yield  complete  self-preserving  solutions  to  the  full  viscous 
equations  of  motion  for  isotropic  turbulence. 

Without  any  loss  of  generality  for  a  complete  self- preserving  isotropic  turbulence  we 
may  set  L  =  X  and  then  substitute  (12)  -  (13)  into  (9)  -  (10),  respectively,  to  get 

Sk  —  —^"{0)  =  constant 

G  =  /‘''(O)  =  constant 

where  a  prime  denotes  a  derivative  with  respect  to  g  =  r/X.  Consequently, 

=  G  =  a  (18) 


where  the  notation  (•)o  denotes  the  initial  value.  The  substitution  of  (18)  into  (6)  and 
(8)  then  yields  the  transport  equations 

K  =  ~e  (19) 

for  complete  self-preserving  isotropic  turbulence.  This  is  a  closed  system  for  the  deter¬ 
mination  of  K  and  e  once  Ko,(-o,Sko  Go  are  provided.  To  simplify  the  subsequent 
presentation,  the  quantity  Go  -  which  is  the  coefficient  of  the  destruction  of  enstrophy 
term  in  (20)  -  will  henceforth  be  referred  to  as  the  “palinstrophy  coefficient”  following 
the  terminology  used  by  Lesieur  (1990).  Accordingly,  the  assumption  of  complete  self- 
preservation  is  seen  to  lead  to  closure  in  the  following  sense;  if  initial  conditions  for  the 
skewness  and  the  palinstrophy  coefficient  are  provided  —  in  addition  to  initial  conditions 
for  K  and  t  —  then  the  energy  decay  can  be  calculated  explicitly  for  all  later  times. 

For  complete  self-preserving  isotropic  turbulence,  the  Karman-Howarth  equation  (5) 
takes  the  form 


5 


or,  equivalently, 


10/  +  27?  —  h?  j-  +  V 


i)  =  R. 


after  replacing  A  using  (11),  (19)  and  (20)  with  R\  =  [u^yi'^X/u  =  ^20/3  ■  Equation 

(22)  will  have  a  solution  if  R\  =  constant  as  first  noticed  by  Dryden  (1943);  this  is  a 
K  ~  t~^  power  law  decay.  However,  (22)  also  has  solutions  where  Rx  =  R\{t)  when 
separability  is  invoked.  The  separability  condition  implies  that  each  side  of  (22)  is  equal 
to  zero  individually,  yielding  differential  equations  from  which  explicit  solutions  for  / 
and  k  may  be  determined  depending  on  the  choice  of  Sko  ^-nd  Go-  These  solutions  were 
first  discovered  by  Sedov  (1944)  and  later  compared  with  experimental  data  by  Korneyev 
and  Sedov  (1976).  The  particular  case  for  which  Go  =  3  so  that  /  is  Gaussian  —  which 
formally  corresponds  to  the  final  period  of  decay  —  was  considered  in  detail  by  Bernard 
(1985).  We  will  briefly  examine  (22)  later  to  establish  the  consistency  of  the  present 
results  with  those  of  previous  studies.  However,  our  analyses  will  be  based  on  a  fixed 
point  analysis  and  direct  numerical  integration  of  (19)  -  (20).  This  will  allow  us  to 
consider  small  departures  from  a  self-preserving  state  as  will  be  demonstrated  later. 


3.  FIXED  POINT  ANALYSIS  AND  NUMERICAL  RESULTS 


In  order  to  carry  out  a  fixed  point  analysis  of  (19)  and  (20),  we  will  combine  these 
equations  into  a  single  transport  equation  for  the  turbulence  Reynolds  number  Rt-  Since, 

Rt  = - ye  (23) 

vt 


it  follows  that 


ij,  =  +  IgA 

u  3\/15  ''  7/  15  1/ 


If  the  transformed  dimensionless  time  r  —  defined  by  the  relation  dr  =  {e/K)dt  —  is 
introduced  into  (24),  we  obtain  the  equation 

The  fixed  points  of  (25)  are  obtained  by  setting  clRt/dr  =  0  which  yields  the  equation; 
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where  (•)co  denotes  the  equilibrium  value  in  the  limit  as  r  — ^  oo.  Equation  (26)  has  the 
solutions: 


co=0 

(27) 

\  2 

(28) 

for  -^Go  <  2,  and 


for  y^Go  >  2.  It  is  a  simple  matter  to  show  that  the  fixed  points  (27)  -  (28)  are  stable 
nodes  that  attract  all  initial  conditions  and  Cq.  It  is  also  evident  from  (28)  that  in  order 
to  have  an  equilibrium  high-Reynolds-number  isotropic  flow  field^  under  self-preserving 
conditions  it  is  necessary  that  Go  ~  \/Rtoo- 

The  fixed  point  Rt^  =  0  is  associated  with  asymptotic  solutions  of  K  and  e  that 
satisfy  the  differential  equations 

K  =  -c  (29) 

(30) 

Equations  (29)  -  (30)  have  the  exact  solution 

where  a  =  IJ^^Go  —  1)  >  1.  For  large  eddy  turnover  times  Cotj Ko  ^  1,  we  have 


K  ~  r“  (33) 

Rt  ~  r“+‘  (34) 

(and,  hence.  K  goes  to  zero  faster  than  Rt),  so  it  follows  that  (31)  -  (32)  can  only  be 
asymptotic  solutions  in  the  limit  as  K  — >  0  and  t  -y  oo.  Therefore,  within  the  framework 
of  complete  self-preservation,  an  asymptotic  power  law  decay  where  K  ~  i~°'  (with 
a  >  1)  is  associated  with  the  fixed  point  Rt^  =  0  and  is  only  formally  consistent  with 
the  final  period  of  decay.  This  is  borne  out  in  the  subsequent  comjjutations. 

During  the  final  period  of  decay  there  is  considerable  evidence  (Batchelor  and  Townsend 
1948a)  indicating  that  /(?/)  =  exp  (-7/^/2)  (i.e.,  a  Gaussian)  in  which  case  (iO)  implies 

^By  a  high  Reynolds  number  isotropic  turbulence  we  mean  the  case  where  R,  >  1;  for  a  low  Reynolds 
number  isotropic  turbulence,  Rt  =  0(1). 
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that  Go  =  3  and,  consequently,  that  a  =  5/2.  The  same  result  is  also  reached  by 
assuming  constancy  of  the  Lmtsianskii  integral 


.  y  oo 

'  /  f{r,t)dr  =  constant 
Jo 


during  decay,  which  appears  to  be  an  acceptable  ;^«s"mption  for  the  final  period.  In 
particular,  from  (12^  it  follows  that  (35)  is  equivalent  to  (c.f.  Hinzc  1975) 


fOO 

Jo 


constant 


so  t'^at  (11)  and  (3C)  imply  that 


=  constant. 


When  (37)  is  combined  with  (31)  and  (32)  it  follows  that  q  =  5/2  —  the  celebrated 
Batchelor  (1948)  result. 

Now  we  will  show  that  the  non-zero  fixed  point  (28)  is  consistent  with  high-Reynolds- 
number  isotropic  turbuler.ee.  The  substitution  of  (28)  into  (19)  -  (20)  yields  the  equations 

k  =  -e  (38) 


e  =  -2- 


K-^Ko  1  + 


tot\  * 


which  have  the  exact  solution 


(namely,  a  power  law  decay  for  the  turbulent  kinetic  energy).  Since  it  can  be  shown 
that  (see  Batchelor  and  Townsend  1948b) 

(42) 

I  e/ h 

it  follows  th^t  G  is  a  ratio  of  turbulent  to  di.s.sipative  time  scales.  It  is  a  simple  matter 
to  show  that  Gee  Hinze  1975) 

-  fTC-  ,2  L'f  .  4\J  -  ^ 

^  ip  J  Q  fv  ^  ^  i  j  cC  K 


where  i)  is  the  energy  spectrum  and  k  is  the  wavenumber.  Imr  complete  self- 
preservation,  E{n^  l)  scales  with  the  Taylor  microscale,  from  which  it  follows  that 


and,  hence,  that 


9u;t 

dx^  dij  ^ 

LUfccet  K 


G  ~  constant 


(4^; 

(45) 


as  shown  earlier.  However,  (43)  is  a  correlation  dominated  by  the  high  wavenumbers 
(i.e.,  small  scales)  and  it  would  tbcredore  seem  more  reasona'  ’  hat  £’(«,/)  should  scale 
with  the  Kolmogorov  length  scale  Ik  =  if  this  is  the  case,  then 


,  .  du>t 

dxj  dxj  ^ 

ujkiCh  Ik 

and 

(47) 

Ik  t 

Equations  i^45)  and  (41)  appear  to  be  contradictory,  however,  they  are  not,  since  in  a 
K  ~  G'-  power  law  decay 

Rt  =  constant.  (48) 

Furthermore,  /?«  =  constant  resolves  the  imbalance  between  the  two  terms  on  the  right- 
hand  side  of  the  dissipation  rate  equation  (20)  since  the  first  teriri  (i.e.,  the  vortex  stretch¬ 
ing)  is  i.otially  of  0{\/Tii)  wh''"'  the  second  term  (i.e.,  the  viscous  destruction)  is  of  0(1). 
This  leads  us  to  the  following  physical  interpretation:  A  K  ~  power  law  decay  is 
the  asymptotic  state  toward  which  a  self-prcserviny  isotropic  turbulence  is  driven  at  high 
Reynolds  numbers  in  order  to  resolve  the  fundamental  imbalance  between  vortex  stretching 
.  nd  viscous  diffusion.  Li  the  process  of  resolving  this  imbalance,  compatibility  with  Kol¬ 
mogorov  scah.ng  is  achieved  for  the  smalt  scale  correlations.  Since  this  consistency  with 
Kolmogorov  scaling  -  which,  on  physical  grounds  should  be  satisfied  at  high  Reynolds 
iiumbers  —  is  achieved  when  ^Go  ~  >  2.  it  is  clear  that  the  non-zero  fixed  point 

(28)  is  the  asymptotic  solution  for  high-Reynolds-uumber  self-preserving  flows. 

Wc  will  now  examine  numerical  solutions  of  (19)  -  (20)  for  complete  self-preservation. 
With  the  possible  exception  of  the  recent  work  of  Bernard  (1985).  we  have  not  seen 
detailed  numerical  re:nilts  published  on  th*'  d^cay  of  K  and  e  in  self-preserving  isotropic 
turbulence.  An  examination  of  these  results  ,viil  ampnfy  the  points  discussed  in  this 


•.ection  and  will  shed  some  interesting  new  light  .m  how  the  self-preservation  assumption 
compares  with  experiments.  In  Figures  1(a)  -  (c),  the  decay  of  the  turbulent  kinetic 
energy  is  shown  for  three  initial  turbulence  Reynolds  numbers  {Rto  —  1000,5000  and 
10,000)  and  three  different  initial  conditions  for  G  (i.e.,  =  1.92,5.0,  and  8.0  where 

Cej  =  ^Go)-  Fot  these  calculations,  as  w^ll  as  the  ones  to  follow,  K'  =  KjKoG'  = 
tJlKo  and  the  skewness 

=  0.5  (49) 

which  is  in  clc-c  proximity  to  the  values  obtained  from  phvs’cai  experiments  in  this  range 
of  Reynolds  numbers.  From  these  figures  it  is  clear  that  the  self-preserving  solution  has 
an  initial  transient  where  the  kinetic  energy  is  fairly  flat;  then  the  kinetic  energy  begins 
to  asymptote  from  above  to  a  power  law  decay  as  evidenced  by  a  straight  line  on  these 
logarithmic  p'ots.  Two  conclusions  can  be  drawn  from  these  results.  First,  for  0^  <  2 
and  Rt^  '>  1,  the  kinetic  energy  does  not  asymptote  to  a  power  law  decay  until 
after  an  extremely  large  number  of  eddy  turnover  times  by  which  time  the  turbulence 
has  decayed  to  a  tiny  fraction  of  its  initial  intensity.  Second,  for  >  2,  the  kinetic 
energy  asymptotes  a  power  law  decay  within  a  few  eddy  turnover  times;  however 
unless  C.j  ~  t/ Rt^  (consistent  with  Kolmogorov  scaling),  the  turbulence  intensity  will 
drop  precipitously  before  this  asymptotic  state  is  achieved. 

To  further  illustrate  these  points,  the  computed  turbulent  kinetic  energy  i^’  compared 
wdth  its  corresponding  asymptotic  power  law  decay  for  increasing  values  of  Figures 

2(a;  -  (c)  which  are  for  initial  turbulence  Reynolds  numbers  of  1000,  10,000  and  100,000, 
respectively.  It  is  clear  from  these  figures  that  for  Cc2  <  2,  the  lurbulent  kinetic  energy 
does  not  reach  its  asymptotic  power  law  dr'cay  even  after  100  eddy  turnover  times!  How¬ 
ever,  for  CVj  sufficiently  larger  than  2,  the  turbulent  kinetic  energy  a.symptotes  to  a 
power  law  decay  within  a  few  eddy  turnover  times.  This  can  be  seen  even  more  vi-  idly  in 
Figure  3  where  q*  =  (i(log /v  *)/ d(log  T)  is  plotted  as  a  function  of  log<‘  for  Rt^  =  1000. 
If  there  is  an  .symptotic  power  law  decay,  this  derivative  will  asymptote  to  the  exponent 
of  the  decay  law.  It  is  clear  that  for  0,2  =  ^ri  exponent  of  1  is  approached  quickly; 
however,  for  Chj  =  F8C  and  1.92  (initial  conditions  which  ultimately  yield  a  power  law 
decay  with  an  exponent  of  approximately  1.2  and  1.1,  respectively)  an  asymptotic  state 
is  not  achieved  even  after  100  edi'y  turnover  times.  Furthermore,  for  <  2  and  large 
initial  turbulence  Reynolds  numbers  Rt^  I,  there  is  a  precipitous  drop  in  the  turbulent 
kinetic  energy  before  a  power  law  decay  is  achiewed;  this  is  due  to  the  early  transient 
when  vortex  stretching  causes  a  considerable  rise  in  tiie  dissipation  (see  Figures  4  (a)  - 
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(b)). 

Since  the  self-preserving  solutions  for  <  2  only  asymptote  to  a  power  law  decay 
in  the  limit  as  K  — >  0  and  t  — ♦  oo,  it  is  reasonable  to  associate  them  exclusively  with 
the  final  period  of  decay.  Experiments  tend  to  indicate  that  the  final  period  of  decay  is 
entered  for  Rt  <  \  wherein  the  exponent  of  the  decay  is  approximately  2.5  (c.f.  Hinze 
1975).  As  noted  earlier,  this  decay  law  is  obtained  asymptotically  for  self-preserving 
isotropic  turbulence  if  C,.,  =  1.4  —  a  result  obtained  by  invoking  Loitsianskii’s  invariant. 
In  Figure  5,  the  decay  of  the  turbulent  kinetic  energy  when  =  1.4  is  shown  for  the 
initial  turbulence  Reynolds  numbers  =  0.25  and  1.0.  It  is  clear  from  these  results 
that  for  Ri^  <  1  the  solutions  begin  to  asymptote  to  the  power  law  decay  (31)  with  an 
exponent  a  =  2.5.  Since  must  equal  1.4  during  the  final  period  of  decay  —  and,  since 

is  a  constant  for  self-preserving  isotropic  turbulence  —  it  follows  that  the  entire  decay 
process  from  high  Reynolds  number  initial  conditions  to  the  final  period  of  decay  cannot 
be  described  within  the  framework  of  complete  self-preservation.  This  conclusion  results 
from  the  fact  that  the  only  consistent  self-preserving  solution  at  high  initial  Reynolds 
numbers  yields  a  K  ~  power  law'  decay  wherein  Rt  asymptotes  to  a  constant  —  a 
state  of  affairs  that  precludes  the  description  of  the  later  stages  of  decay.  Furthermore, 
the  value  of  =  1-4,  which  describes  the  final  period  of  decay,  yields  unphysical  results 
for  the  early  stages  of  a  high  Reynolds  number  isotropic  turbulence  (i.e.,  it  predicts  an 
early  time  transient  where  there  is  a  precipitous  drop  in  the  turbulent  kinetic  energy; 
see  Figure  6).  In  order  to  describe  the  entire  decay  process  of  a  high  Reynolds  number 
isotropic  turbulence,  G  as  well  as  5k  must  vary  with  time  —  a  possibility  that  is  precluded 
by  the  assumption  of  self-preservation  which  renders  them  to  be  constant. 

4.  COMPARISONS  WITH  ALTERNATIVE  THEORETICAL  ANALYSES 
AND  EXPERIMENTS 

The  results  derived  in  the  pre”ious  section  are  consistent  with  those  of  Batchelor 
(1948)  for  low  turbulence  Reynolds  numbers;  however,  our  high  Reynolds  number  asymp¬ 
totic  solution  yields  K  ~  t"'  whereas  in  Batchelor’s  solution  K  ~  reason 

for  this  difference  is  simple:  as  alluded  to  earlier,  Batchelor  also  found  the  K  ~  t~^ 
solution  but  dismissed  it  as  a  viable  result  since  Loitsianskii’s  integral  was  not  a  dy¬ 
namic  invariant.  Interestingly  enough,  an  earlier  experimental  study  by  Batchelor  and 
Townsend  (1948b)  yielded  results  that  were  far  more  suggestive  of  a  K  -  f*  rather 
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than  a  K  ~  power  law  decay.  Despite  the  fact  that  Batchelor  (1948)  states  that 

the  K  ~  decay  law  is  a  complete  self-preserving  solution,  in  reality  it  is  only  a 

partial  self-preserving  solution  corresponding  to  the  inviscid  Karman-Howarth  equation 
(see  Hinze  1975).  It  is  our  view  that  since  the  K  ~  asymptotic  decay  law  is  a  formal 
solution  to  the  full  Karman-Howarth  equation,  it  should  not  be  dismissed  unless  it  is  in 
incontrovertible  contradiction  of  experiments  or  other  exact  theoretical  results. 

As  mentioned  earlier,  Dryden  (1943)  postulated  a  A"  ~  power  law  decay  based 
on  a  direct  analysis  of  the  Karman-Howarth  equation.  He  observed  —  as  is  evident  from 
(22)  —  that  the  Karman-Howarth  equation  will  allow  for  self-similar  solutions  if 


R\  =  constant 


(50) 


which  yields  K  ~  t~^  as  a  direct  consequence  of  (6).  However,  there  are  other  temporally 
varying  solutions  to  (22);  complete  self-preservation  only  requires  that  R\  asymptote  to 
a  constant.  Sedov  (1944)  studied  solutions  of  (22)  obtained  by  applying  the  separability 
constraint 


df\ ,  jf  n 


di]  \  dr] )  dr] 


3^0  -  5  )=  0 


=  0 


(51) 


(52) 


which  renders  R\  =  R\{t)  consistent  with  (19)  -  (20).  Solutions  to  (51)  -  (52)  have 
not  been  studied  in  great  depth  subsequent  to  Sedov  (1944)  who  showed  that  /(?;)  = 
M(^,  2^)  where  M  is  the  confluent  hypergeometric  function  and  7  =  |Go  — 

Batchelor  (1948)  expressed  concern  over  the  fact  that  this  solution  leads  to  a  unique 
determination  of  both  /  and  k\  however,  although  he  suspected  that  the  Sedov  solution 
was  unphysical,  he  stated  that  he  was  “not  able  to  find  any  definite  anomalies”.  In 
the  limit  as  R\  — >  0  it  is  clear  that  (51)  is  a  direct  consequence  of  the  Karman-Howarth 
equation.  Consequently,  it  is  not  surprising  that  the  Sedov  solution  for  the  final  period  of 
decay  yields  physically  interesting  solutions  as  recently  demonstrated  by  Bernard  (1985). 
However,  it  will  now  be  shown  definitively  that  the  Sedov  solution  yields  unphysical 
results  at  high  turbulence  Reynolds  numbers.  In  figure  7,  the  results  of  a  numerical 
solution  of  (51)  for  /  are  shown  for  a  variety  of  values  of  Go  ranging  from  3  to  60.  For 
Go  =  3  it  can  be  shown  analytically  that  /  =  exp(— t/^/2)  yielding  an  energy  spectrum 
of  the  form 

^*(«*)  =  exp(-K*V2)  (53) 
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where  k*  =  kX  and  E*  =  Eju^X.  Equation  (53)  is  obtained  from  the  identity 


E’[k’')  =  —I  f{r]){K*r}  sin  k^t]  —  k*  t]^  cos  K'‘r))d7]  (54) 

TT  Jo 

(c.f.  Batchelor  1953).  This  result  —  which  has  £’*(«’)  ~  at  low  wavenumbers  and 
has  £’(«*)  decaying  exponentially  at  high  wavenumbers  —  is  consistent  with  established 
results  on  the  final  period  of  decay  (c.f.  Hinze  1975).  However,  for  sufficiently  large  Go, 
it  is  a  simple  matter  to  show  from  (51)  that 

/(,)  ~  (55) 


for  7/  >•  1.  This  explains  why  /{rj)  is  so  slow  to  asymptote  to  zero  when  Go  >  10  in 
figure  7.  In  fact  for 

-a  >2 


it  follows  from  (54)  and  (55)  that  E’“(k*)  becomes  singular.  From  (28)  it  can  then  be 
concluded  that  the  Sedov  solution  will  yield  a  singular  energy  spectrum  when 


Rtoo  >  0 

It  is  thus  clear  that  the  Sedov  solution  yields  untenable  results,  at  high  Reynolds  numbers, 
for  the  double  and  triple  two-point  longitudinal  velocity  correlations. 

The  major  deficiency  with  the  approaches  of  Dryden  and  Batchelor  —  as  well  as  that 
of  Sedov  —  lies  in  the  use  of  the  self-similar  Karman-Howarth  equation  (22)  which  does 
not  allow  for  the  treatment  of  small  departures  from  a  state  of  complete  self-preservation. 
Such  small  departures  can  be  characterized  by  the  perturbations 


SK  =  SKo  +  SSKit)  (56) 

G  =  Go  +  SG{t)  (57) 

where  1|  6Sk  ||  /Sko  <  T  II  1|  /Go  <  1  and  6SKit),^G{t)  ^  0  as  t  ^  oo.  The 

substitution  of  (56)  -  (57)  into  (6)  and  (8)  yields  the  governing  equations  for  small 

departures  from  a  state  of  self-preservation.  If  we  denote  SK  and  Se  as  the  departures 
from  the  self-preserving  solutions  K  and  e  obtained  from  (19)  -  (20),  it  follows  that  for 
the  perturbations  (56)  -  (57)  we  will  have 


II  II 
II  i<  II 


<  1, 


liM 

lull 


« 1 
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due  to  the  fact  that  (19)  -  (20)  have  fixed  points  that  are  stable  nodes  (c.f.  Guckenheimer 
and  Holmes  1986).  Consequently,  (19)  -  (20)  will  yield  an  excellent  approximation  for 
isotropic  decay  when  there  are  extremely  small  departures  from  a  state  of  complete  self- 
preservation.  In  contrast  to  this  nice  behavior,  the  Karman-Howarth  equation  becomes 
indeterminate  when  subjected  to  infinitesimal  perturbations  from  a  self-preserving  state. 
Hence,  it  appears  that  the  one-point  equations  (19)  -  (20)  form  a  broader  basis  for  the 
analysis  of  the  energy  decay  of  self- preserving  isotropic  decay  than  does  (22). 

The  general  solution  to  the  complete  self-preserving  isotropic  decay  equations  (19)  - 
(20)  at  high  Reynolds  numbers  is  shown  schematically  in  figure  8  (this  is  for  the  physically 
significant  case  where  >  Rt^  so  that  the  turbulence  Reynolds  number  decays).  There 
is  an  early  time  transient  (region  AB)  where  the  turbulent  kinetic  energy  is  flat;  it  is 
eventually  followed  by  the  asymptotic  region  CD  where  K  ~  .  These  two  regions  are 

connected  by  the  overlap  region  BC.  The  initial  transient  AB  evolves  on  the  Kolmogorov 
time  scale  during  which  time  there  is  a  precipitous  drop  in  the  turbulence  Reynolds 
number  (see  figure  9).  On  the  other  hand,  the  overlap  region  BC  evolves  on  the  turbu¬ 
lence  time  scale  Kjc,  in  this  region  the  turbulence  Reynolds  number  Rt  becomes  close  to 
approaching  it  asymptotically  from  above.  As  a  direct  consequence  of  the  perturba¬ 
tion  analysis  discussed  above,  the  overlap  region  BC  can  be  set  into  strong  approximate 
agreement  with  the  asymptotic  approach  to  a  state  of  complete  self-preservation.  These 
results  have  a  direct  bearing  on  how  the  complete  self-preserving  solution  compares  with 
physical  experiments  as  we  will  soon  see. 

It  is  widely  believed  that  a  K  ~  asymptotic  decay  law  is  in  violation  of  experi¬ 
mental  data  for  isotropic  turbulence.  This  experimental  data  (see  Uberoi  1963,  Kistler 
and  Vrebalovich  1966,  Comte-Bellot  and  Corrsin  1966,  1971,  and  Warhaft  and  Lumley 
1978)  has  yielded  power  law  decays  with  exponents  varying  from  1  to  1.4  with  a  mean  of 
approximately  1.25.  However,  great  caution  must  be  taken  in  using  this  data  to  dismiss 
the  possibility  of  a  K  ~  asymptotic  power  law  decay  at  high  Reynolds  numbers  since 
most  of  this  data  is  for  a  limited  number  of  eddy  turnover  times  (typically  for  totjKo  <  4). 
If  the  self-preserving  solution  is  examined  for  this  same  limited  number  of  eddy  turnover 
times  it  follows  that  the  resulting  solution  can  be  fitted  to  an  excellent  degree  of  approx¬ 
imation  by  a  power  law  decay  with  exponents  in  the  range  of  1  to  1.4  depending  on  the 
initial  conditions;  the  lower  the  Reynolds  number,  the  longer  the  solution  takes  to  reach 
an  asymptotic  state  and  the  larger  the  exponent  is  during  the  early  stages  of  decay  (see 
Figures  10  (a)  -  (c)).  Consequently,  if  one  examined  in  isolation  the  self-pre.serving  solu- 
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tions  for  the  first  few  eddy  turnover  times  (with  the  short  early  time  transient  omitted), 
one  could  erroneously  conclude  that  there  was  an  asymptotic  power  law  decay  with  an 
exponent  in  the  range  of  1  to  1.4  depending  on  the  initial  conditions;  in  reality,  all  of 
these  solutions  are  asymptoting  to  a  power  law  decay.  The  solutions  shown  in  Figures 
10  (a)  -  (c)  correspond  to  the  overlap  gion  BC  shown  in  Figure  8  and,  hence,  can  be 
associated  with  the  asymptotic  approach  to  a  state  of  complete  self-preservation.  .An 
argument  has  been  raised  recently  by  Walker  and  Corrsin  (1985)  and  Walker  (1986)  that 
the  physical  experiments  may  not  go  far  enough  to  see  a  power  law  decay.  Unless 
the  initial  turbulence  Reynolds  number  is  extremely  large,  an  asymptotic  state  may  not 
be  achieved  in  the  first  few  eddy  turnover  times.  In  this  regard  it  is  interesting  to  note 
that  the  only  extremely  high-Reynolds-number  experiment  (i.e.,  Kistler  and  Vrebalovich 
1966)  and  large-eddy  turnover  time  experiment  (Walker  1986)  did  measure  a  K  ~ 
asymptotic  power  law  decay.  Consequently,  existing  experiments  cannot  rule  out  the 
possibility  of  a  K  ~  asymptotic  power  law  decay  at  high  Reynolds  numbers  and 
do  not  warrant  the  dismissal  of  the  complete  self-preserving  solution  discussed  herein. 
Furthermore,  Rosen  (1981)  recently  derived  a  asymptotic  power  law  decay  based  on 
alternative  ideas  from  statistical  mechanics. 

Now,  we  will  address  the  interesting  controversy  generated  recently  by  George  (1987, 
1989).  He  claimed  to  find  complete  self-preserving  solutions  —  with  the  Taylor  microscale 
as  the  similarity  length  scale  —  that  exist  for  all  turbulence  Reynolds  numbers.  These 
solutions  were  characterized  by  an  asymptotic  power  law  decay  where  the  exponent  is 
determined  by  the  initial  conditions.  George  arrived  at  this  alternative  self-preserving 
solution  by  relaxing  the  classical  similarity  constraint  (13).  He  argued  that  the  normal¬ 
ization  of  the  two-point  triple  velocity  correlation 


T{r,t)  —  u2(x,  t)u{x  -f  r,  t) 


(58) 


by  to  form  k{r,  t)  is  arbitrary  since  its  one-point  contraction  7’(0,  t)  =  is  zero. 

(This  stands  in  contrast  to  the  formulation  of  /(r,  t)  which  is  obtained  by  normalizing  the 
two-point  double  velocity  correlation  with  its  one-point  contraction  li^).  Consequently, 
George  argued  that  constraint  (13)  should  be  replaced  with  the  alternative  constraint 


Tir,t)  ^  -  /r \ 
w{t)  \\)' 


(59) 


where  'w{t)  is  a  suitable  weighting  function.  Then  —  from  the  definition  of  Sj^  in  equation 


(9)  and  the  Karman-Howarth  equation  (22)  —  instead  of  the  constraints 

R)^  =  constant^  Sk  =  constant  (60) 

which  renders  a  t~^  power  law  decay,  we  get  the  constraints 

R\Sk  =  constant,  R\wK~^^^  =  constant  (61) 


which  allows  for  the  possibility  of  an  alternative  decay  law  within  the  general  framework 
of  complete  self-preservation.  From  (61)  it  follows  that 

Sk  oc  (62) 


and  that 

w  «  (63) 

which  shows,  incidentally,  that  w  cannot  be  chosen  arbitrarily.  Since  the  proportionality 
constant  in  (62)  must  be  Sna^ti^  ~  and  since  G  still  remains  a  constant  Go  during  the 
decay  —  this  alternative  self-preservation  leads  to  the  decay  equations 

=  -e  (64) 


—  77  ~  77^077 


(65) 


A  7  ^2 

7  d1/2  ^ 

3x/T5  K 

instead  of  (19)  -  (20).  Equations  (64)  -  (65)  yield  the  closed  form  solution  for  the  energy 
decay 


-0 


(66) 


where 

This  is  indeed  a  power  law  decay  with  an  exponent  that  depends  on  the  initial  conditions 
as  claimed  by  George  (1987,  1989). 

Although  we  do  feel  that  George  raised  some  interesting  issues  —  and  made  a  seem¬ 
ingly  valid  point  concerning  the  arbitrariness  of  the  normalization  of  the  two-point  triple 
velocity  correlation  —  it  is  our  opinion  that  there  is  a  problem  with  this  alternative 
self-preserving  solution.  By  virtue  of  (61),  Geoige’s  alternative  scaling  renders  a  zero 
skewness  in  the  limit  of  infinite  turbulence  Reynolds  numbers.  Experiments  indicate 
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that  the  skewness  is  an  absolute  constant  of  order  one  at  extremely  high  Reynolds  num¬ 
bers  (see  Van  Atta  and  Antonia  1980).  Furthermore,  a  direct  physical  consequence  of  a 
vanishing  skewness  Sk  is  the  occurrence  of  zero  transfer  and,  hence,  no  energy  cascade. 
This  can  be  seen  easily  from  equation  (59)  which  can  be  re-written  as 

T{r,t)  oc  (^0 

which  implies  that  the  transfer  T(r,t)  0  as  Rt  oo.  Vanishing  transfer  in  the  limit 
of  infinite  Reynolds  numbers  is  an  untenable  physical  result  and  hence  this  alternative 
self-preserving  solution  is  not  acceptable  for  high-Reynolds-number  isotropic  turbulence. 

Finally,  a  few  comments  are  in  order  concerning  the  implications  of  these  results 
for  turbulence  modeling.  In  the  commonly  used  turbulence  models,  the  dissipation  rate 
equation  is  modeled  as 

^  =  (68) 

for  isotropic  decay,  where  is  a  constant  (c.f.  Launder  and  Spalding  1974,  Speziale 
1991).  Equation  (68)  is  derived  by  invoking  Kolmogorov  scaling  for  G  which  requires 
that 

G  =  C’iV^-hC2  (69) 

where  Ci  and  C2  are  constants;  an  equilibrium  hypothesis  is  then  made  by  which  C\  = 
SO  that  the  leading  order  part  of  the  destruction  of  dissipation  term  annihilates 
the  vortex  stretching  term  in  (8)  yielding  (68).  In  contrast  to  (68),  the  complete  self¬ 
preserving  solution  has  a  dissipation  rate  transport  equation  of  the  general  mathematical 


where  and  are  constants.  Equation  (70)  can  be  also  derived  based  on  Kolmogorov 

scaling  (69)  when  departures  from  equilibrium  are  allowed  wherein  C\  The 

addition  of  the  unbalanced  vortex  stretching  term  in  (70)  allows  for  a  better  treatment 

of  departures  from  equilibrium  in  several  ways.  First,  as  shown  earlier  in  Figures  10  (a)  - 

(c),  the  self-preservation  model  allows  for  the  description  of  the  initial  stages  of  isotropic 

decay  where  the  exponent  of  the  decay  law  can  vary  mildly  with  the  initial  conditions  — 

a  feature  observed  in  physical  experiments.  On  the  other  hand,  the  more  commonly  used 

1 

model  (68)  predicts  a  universal  decay  law  where  K  ~  i  for  all  and  all  t  >  0. 

Second,  the  self-preserving  solution  can  accommodate  the  limit  of  zero  viscosity.  In  this 
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limit,  it  is  a  simple  matter  to  show  that  (20)  predicts  a  finite  time  enstrophy  blow-up  at 
the  critical  time 

.  _  6^15 

where  is  the  initial  enstrophy.  In  Figure  11(a),  the  time  evolution  of  the  enstrophy 
corresponding  to  the  complete  self-preserving  solution  is  shown  for  a  variety  of  increasing 
Rt„]  it  is  clear  that  an  enstrophy  crisis  is  predicted  for  Rt„  >■  1  which  eventually  leads 
to  a  finite  time  enstrophy  blow-up  in  the  limit  as  i/  — >  0.  These  results  are  in  excellent 
qualitative  agreement  with  results  obtained  from  EDQNM  as  illustrated  in  Figure  11(b) 
taken  from  Lesieur  (1990).  While  the  issue  of  a  finite  time  enstrophy  blow-up  is  still 
being  debated  by  the  turbulence  community  (c.f.  Pumir  and  Siggia  1990),  one  thing  is 
clear:  the  enstrophy  grows  dramatically  when  i/  =  0.  In  contrast  to  the  results  shown  in 
Figures  11  (a)  -  (b),  the  commonly  used  dissipation  rate  model  (68)  erroneously  predicts 
that  the  enstrophy  is  conserved  in  the  inviscid  limit,  i.e.  that 

=  constant  (72) 

when  1/  =  0.  It  thus  appears  that  the  complete  self-preserving  solution  allows  for  a 
better  treatment  of  non-equilibrium  isotropic  turbulence  that  could  be  of  future  use  in 
the  development  of  improved  turbulence  models. 

5.  CONCLUSIONS 

The  energy  decay  for  complete  self- preserving  isotropic  turbulence  has  been  re-examined 
from  a  basic  theoretical  and  computational  standpoint.  Several  interesting  conclusions 
can  be  drawn  from  these  results: 

(1)  The  nonlinear  differential  equations  for  the  energy  decay  have  two  fixed  points  — 

Ri^  =  0  and  Rt^  =  ~  The  former  fixed  point  is  only  achieved  in  the 

limit  as  t  — >  oo  and  hence  is  associated  with  the  final  period  of  decay.  Consistent  with 
the  Batchelor  (1948)  result,  a  K  ~  t~^^^  power  law  decay  is  obtained  when  Loitsianskii’s 
invariant  or  the  Gaussianity  of  /(?/)  is  invoked. 

(2)  The  non-zero  fixed  point  Rt^  =  ^{^Go  —  2)^/5^'^  is  approached  within  a  few 
eddy  turnover  times  and  gives  rise  to  a  ~  t~^  asymptotic  power  law  decay.  It  is  the 
high-Reynolds-number  asymptotic  solution  for  a  complete  self-preserving  isotropic  tur¬ 
bulence.  This  solution  appears  to  have  been  prematurely  dismissed  by  Batchelor  (1948) 
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purely  on  the  grounds  that  Loitsianskii’s  integral  was  not  an  invariant  —  a  constraint 
which  wa^  later  learned  to  be  violated  in  isotropic  turbulence  when  Rt^  ^  1. 

(3)  The  structure  of  the  high  Reynolds-number  self-preserving  solution  during  the 
first  few  eddy  turnover  times  was  examined  in  detail.  By  a  perturbation  analysis,  it  was 
argued  that  these  solutions  can  serve  as  an  approximation  for  the  asymptotic  approach 
to  a  state  of  complete  self-preservation.  It  was  found  that,  depending  on  the  initial 
conditions,  the  early  time  solutions  could  be  fitted  with  a  power  law  decay  which  has  an 
exponent  varying  from  1.0  to  1.4  —  a  range  of  values  that  is  compatible  with  existing 
experimental  data.  Consequently,  existing  experiments  cannot  rule  out  the  possibility 
of  a  complete  self-preserving  solution  with  a  K  ~  asymptotic  power  law  decay  at 
high  Reynolds  numbers.  In  fact,  this  asymptotic  decay  law  was  found  recently  in  the 
experiments  of  Walker  (1986). 

(4)  Since  the  assumption  of  complete  self-preservation  requires  that  G  be  constant 

—  and  since  for  high-Reynolds-number  isotropic  turbulence  —G  >  2,  whereas  for  low- 
Reynolds- number  isotropic  turbulence  <  2  —  ii  is  clear  that  the  entire  process  of 
isotropic  decay  from  high  Reynolds  number  initial  conditions  to  the  final  period  of  decay 
cannot  be  described  by  the  theory.  The  alternative  conclusion  drawn  by  George  (1987, 
1989)  arose  from  the  use  of  a  new  scaling  for  the  two-point  triple  velocity  correlation  that 
rendered  a  time  dependent  skewness  where  Sk  ~  However,  this  type  of  asymptotic 

behavior  for  the  skewness  wherein  Sk  — >  0  as  >  oo  is  physically  inconsistent  since  it 
yields  zero  energy  transfer  in  the  infinite  Reynolds  number  limit. 

Within  the  framework  of  self-preservation,  the  physical  origin  of  a  power  law 

decay  becomes  clear:  it  is  the  asymptotic  state  toward  which  a  high-Reynolds  number 
isotropic  turbulence  is  driven  in  order  to  resolve  an  0{R/  )  imbalance  between  vortex 
stretching  and  viscous  diffusion.  The  resolution  of  this  imbalance  also  yields  compatibility 
with  Kolmogorov  scaling.  Results  were  presented  which  indicate  that  the  complete  self¬ 
preserving  solution  yields  a  better  description  of  non-equilibrium  isotropic  turbulence 
than  the  commonly  used  turbulence  models.  It  is  also  interesting  to  note  that  when 
the  self-preserving  assumption  is  extended  to  homogeneous  shear  flow,  a  production- 
equals-dissipation  equilibrium  can  occur  —  preceded  by  a  transient  where  K  and  e  grow 
exponentially  —  as  recently  shown  by  Bernard  and  Speziale  (1990).  It  thus  appears  that 
the  theory  of  self-preservation  in  homogeneous  turbulence  has  many  interesting  features 
that  have  not  yet  been  fully  understood  and  are  worthy  of  further  study. 
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Figure  1.  Decay  of  turbulent  kinetic  energy  in  complete  self-preserving  isotropic  turbu¬ 
lence  for  initial  turbulence  Reynolds  numbers  =  1000,5000  and  10,000:  (a)  — 

1.92,  (b)  Cc,  -  5,  and  (c)  =  8. 


P'igure  1.  Decay  of  turbulent  kinetic  energy  in  C'-mplete  self-preserving  isotiopir  turbu¬ 
lence  for  initial  turbulence  Reynolds  numbers  Rt^  =  1000,5000  and  10,000:  (a)  C,,  = 
1.92,  (b)  C,,  =  5,  and  (c)  C,,  =  ?. 


Figure  1.  Decay  of  turbulent  kinetic  energy  in  complete  self-preserving  isotropic  turbu¬ 
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Figure  2.  Decay  of  turbulent  kinetic  energy  for  a  variety  of  initial  conditions  on  G 

(Ctj  =  jiGo)'-  -  self-preserving  solution; - asymptotic  solution  K  ~  where 

a  =  {C,,  -  l)-‘  for  C,,  <  2  and  a  =  1  for  C,,  >  2.  (a)  Rt„  =  1000,  (b)  Rt„  =  10,000 
and  (c)  Rt^  =  100,000. 
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Figure  2.  Decay  of  turbulent  kinetic  energy  for  a  variety  of  initial  conditions  on  G 

((7^  _  2-G„):  -  self-preserving  solution; - asymptotic  solution  K  ~  where 

a  =  [Cc2  —  1)”^  for  <  2  and  q  =  1  for  >  2.  (a)  Rt„  —  1000,  (b)  Rt„  =  10,000 
and  (c)  =  100,000. 
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Figure  8.  Schematic  of  the  complete  self-preserving  solution  for  the  decay  of  turbulent 
kinetic  energy  at  high  Reynolds  numbers  {Rt^  >  Rt^). 
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Figure  11.  Time  evolution  of  the  enstrophy  at  high  turbulence  Reynolds  numbers:  (a) 
complete  self-preserving  solution  {Ci^  =  2)  and  (b)  schematic  from  Lesieur  (1990)  based 
on  EDQNM. 
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